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Foreword 


Numerous naval applications are directly dependent on knowing the past, 
present, and future sea surface temperature (SST) fields. The Diurnal Ocean 
Surface Layer (DOSL) model at the Fleet Numerical Oceanography Center 
attempts to fill a portion of these requirements by forecasting the 24-hour 
change in global SSTs. This procedure is done by using the wind stress and 
heat flux forcing functions available from the Naval Operational Global 
Atmospheric Prediction System as input to the parameterized DOSL model. 

This report deals with validating the DOSL products using infrared imagery 
obtained in several locations during the night and day to map the diurnal SST 
changes. Particular attention is focused on extreme diurnal warming events, 
since these cases are most commonly associated with surface mixed-layer 
warming that disrupts the surface acoustic duct. Several events are detailed 
in the subtropical Atlantic to depict the capabilities within the DOSL model. 
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Technical Director 


J. B. Tupaz, Captain, USN 
Commanding Officer 




Executive Summary 


The Diurnal Ocean Surface Layer (DOSL) model at the Fleet Numerical 
Oceanography Center forecasts the 24-hour change in global sea surface 
temperatures (SST). Validating the DOSL model is a difficult task due to the 
huge areas involved and the lack of in situ measurements. Therefore, this report 
details the use of satellite infrared multichannel SST imagery to provide day 
and night SSTs that can be directly compared to DOSL products. This water- 
vapor-corrected imagery has the advantages of high thermal sensitivity 
(0.12°C), large synoptic coverage (nearly 3000 km across), and high spatial 
resolution that enables diurnal heating events to be readily located and mapped. 

Several case studies in the subtropical North Atlantic readily show that DOSL 
results during extreme heating periods agree very well with satellite-imagery- 
derived values in terms of the pattern of diurnal warming. The low wind and 
cloud-free conditions necessary for these events to occur lend themselves well 
to observation via infrared integer y. Thus, the normally cloud-limited aspects 
of satellite imagery do not come into play for these particular environmental 
conditions. 

The fact that the DOSL model does well in extreme events is beneficial from 
the standpoint that these cases can be associated with the destruction of the 
surface acoustic duct. This so-called “afternoon effect” happens as the after- 
noon warming of the mixed layer disrupts the sound channel and the 
propagation of acoustic energy. Thus, the DOSL model can be used to help 
estimate the location of where the afternoon effect may occur, as well as 
help analysts interpret artificially induced, diurnal heating SST gradients when 
generating mesoscale front and eddy maps. 
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Diurnal Ocean Surface Layer Model Validation 


I. Introduction 

The U.S. Navy requires a wide range of real-time 
sea surface temperature (SST) analyses (nowcasts) and 
forecasts to meet the increasingly sophisticated needs 
demanded by sensors and weapon systems. The SSTs 
may directly or indirectly impact successful and 
efficient completion of duties ranging from antisub- 
marine warfare (ASW), optimum track ship routing, 
and search and rescue, or serve as vital input to 
atmospheric prediction models that require accurate 
air-sea interface initial conditions. This report focuses 
on the validation of one model aimed at satisfying the 
need to forecast daily SST changes. 

The Diurnal Ocean Surface Layer (DOSL) model 
(Clancy, 1986) is run daily at the Fleet Numerical 
Oceanography Center (FNOC) in Monterey, CA, and 
predicts the daily SST changes due to the impact of 
surface heat fluxes and wind stress. This parameterized 
model requires no SST initial conditions and uses the 
forecast winds and fluxes from the Navy’s global 
atmospheric model to provide the key input 
parameters. It then produces a prediction for the 
diurnal SST heating that is valid globally at the time 
of local maximum heating. 

The DOSL model thus generates a contoured map 
for both hemispheres outlining the regions of large 
diurnal heating. This graphical guidance can be used 
to outline areas where the surface acoustic duct may 
be lost due to solar-heating-induced stratification 
(afternoon effect); to help correctly interpret infrared 
SST or forward-looking infrared images for ocean 
mesoscale feature locations; and to delineate regions 
where day and night satellite SSTs should not be mixed, 
since they represent two totally different environmental 
conditions. The DOSL model can thus provide some 
unique capabilities when properly used in conjunction 
with the right Navy mission. 

Validating a prognostic SST model over large time 
and space scales (DOSL uses a polar stereo grid of 
~ 320 km and is interpolated to a 2.5° x 2.5° grid) 
raises several interesting problems. It is not sufficient 
to match individual point values with in situ measure- 
ments, since the pattern of SST warming is also a key 
output value. Therefore, accurate SST reports within 
regional and sometimes basin scale domains must be 


combined with the attribute of coincident timing 
needed at the minimum and maximum periods of solar 
heating. Multichannel SSTs (MCSST) from polar 
orbiter infrared sensors are the only means presently 
available to attempt this arduous task. 


II. Diurnal Ocean Surface Layer Model 

The DOSL model uses the scaling relationships 
derived by Price et al. (1986) to relate the expected 
afternoon mixed-layer response to the local 
atmospheric forcing. A model similar to that of Price 
et al. (1978) was used to simulate the extensive ocean 
mixed-layer data set obtained from the R/V Flip using 
rapid-profiling conductivity-temperature-depth and 
VMCM current meters. These detailed observations 
and model runs enabled Price et al. (1986) to 
develop and verify scaling relationships that predict 
the afternoon SST warming and transient layer depth 
using surface heat fluxes and surface wind stress values. 
This parameterization was later verified by Bowers 
et al. (1986) and Price et al. (1987) with independent 
data from the northwestern Atlantic. 

These scaling relationships are the basis for the 
DOSL model now functioning in real time with global 
coverage at FNOC (Clancy, 1986). It forecasts the 
change in day-night SST (AT) and diurnal thermocline 
depth (D), valid everywhere in the local afternoon. This 
diurnal SST variation will occur locally from a 
minimum near 0600 to a maximum about 1400 local 
time. The change due to diurnal heating is given by 
Equation 1: 

AT =1.5 ^ 

C 

where g is the acceleration of gravity, a is the expansion 
coefficient for seawater, C is the heat capacity of 
seawater, H is the maximum net downward surface 
heat flux (latent + sensible + infrared + solar) during 
the day, P is 1/2 the time interval over which the net 
surface heat flux is downward, S is the wind stress 
parameter, and F is a radiation extinction function. 
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The parameter S is calculated from 
S = ^ [1 - cos (fP)]° 5 , (2) 

where t is the surface wind stress averaged over the 
period of net downward surface heat flux and / is 
the Coriolis parameter. The radiation function F 
is derived from 

_ . 0.4 (2/-//) exp (-kD) m 


where / is the maximum surface solar heat flux during 
the day and k is the extinction coefficient for solar 
radiation (Clancy, 1986). 

The forcing functions are derived from the Navy 
Operational Global Atmospheric Prediction System 
(NOGAPS) model (Rosmond, 1981). This model 
provides forcing functions on a global, 73 x 144 
spherical grid that yields a resolution of 2.5°. These 
forcing fields are subsequently interpolated to the 
63 X 63 polar stereographic grid on which DOSL runs. 
This latter grid provides resolutions that vary from 
200 km (tropics) to nearly 400 km (polar regions). In 
determining S, the wind-speed-dependent drag 
coefficient of Garrat (1977) is used to calculate surface 
wind stresses from the NOGAPS 19.5-m winds. Also, 
the algorithms of Clancy and Poliak (1983) are used 
to interpolate the surface heat fluxes and wind stresses 
between the 6-hourly values provided by NOGAPS 
(Clancy, 1986). 

Particular attention should be paid to the forcing 
function grid resolution. The subtropical zones with 
the best chances of strong diurnal heating occur 
between 20° and 30° from the equator. NOGAPS has 
a resolution near 275 km and, thus, cannot clearly 
define mesoscale heating events that may occur and 
cannot be verified with high-resolution, satellite-derived 
SSTs. The forcing function spatial definition is, thus, 
one limiting factor in the DOSL accuracy. 

The DOSL model produces a forecast once per day 
at 00Z (Greenwich Mean Time), which is valid 
worldwide in the local afternoon. This means that the 
valid DTG (Date-Time-Group) (86081500 is 00Z 
15 Aug 86) is a function of longitude. The valid time 
at a particular location is 

Time = DTG + 14 hr + - (lhr) 

15° 

for 0 < L < 270° W (4) 

and 

Time = DTG + 8 hr + (L ~ 270) 

15° 

for 270 < L < 360° W , (5) 


where L is the longitude of the location ranging 
from 0° to 360° measured positive westward (W) from 
Greenwich. 

The table illustrates sample valid forecast times at 
different longitudes for the same DTG of 86081500. 


Longitude 

L (°W) 

Valid DTG (Z) 

Local Time 

15°W 

15° 

86081515 

1400 

75’W 
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to 
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86081519 

1400 

150°W 

150° 

86081600 

1400 
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© 

o 

T“ 

CM 

86081604 

1400 

60° E 

300° 

86081510 

1400 

15°E 

345° 

86081513 

1400 


Sample graphical DOSL output (6 August 1987) is 
shown in Figure 1. The large FNOC chart size and 
day/night SST difference contours at 0.2°C intervals 
necessitate exhibiting just the Atlantic and East Pacific 
basins; otherwise, interpreting values would be quite 
difficult. One can readily locate areas where significant 
warming (greater than 1°C) is forecast by the closely 
bunched contours in numerous regions. 

The North Atlantic is particularly active, with greater 
than 1 .5°C warming predicted by the DOSL model for 
the Gulf of Mexico, the Gulf Stream, the Northeast 
Atlantic south of Iceland, and the two local maxima 
at 30°N off Northwest Africa. The number and 
intensity of changes forecast for 6 August 1987 is not 
the norm. A scan through weeks of charts reveals 
that the patterns change dramatically in size, in shape, 
and in value from one day to the next as the synoptic 
weather pattern evolves. 

Figure 2 illustrates the corresponding surface wind 
and atmospheric pressure fields (6 August 1986) 
available from the FNOC marine wind analysis. The 
1024-mb subtropical high centered between 20°-30°N 
stretches across a huge longitudinal range and readily 
dominates the weather and, thus, cloud cover. This 
event is a critical factor in that subsidence within the 
high will tend to deter cloud development, thus opening 
the way for high insolation. Light winds are also the 
rule, as evidenced by the lack of wind barbs within 
the center. (Higher resolution wind fields will be 
covered later.) 

III. Diurnal Warming and 
SST Measurements 

Insolation affects the entire mixed-layer water 
column to various degrees, depending on several 
factors: amount of solar heating, surface wind stress, 
water clarity, and temperature stratification when the 
heating event first begins. The warming is initially more 
abrupt at the sea surface and is thus reflected in the 
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150°W 140° 130° 120° 110° 100° 90° 80° 70° 60° 50° 40° 30° 20° lO°W 0° 10°E 20° 30° 40° 50° 60°E 

Figure I. Sample graphical output from the DOSL model for 6 August 1987. Contours every 0.2°C represent forecast of diurnal change in SST 
valid everywhere at local time of maximum insolation. (Note all chart values should be multiplied by 0. 1.) 
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150°W 140° 130° 120° 110° 100° 90° 80° 70° 60* 50° 40° 30° 20° 10°W C° 10°E 20° 30° 40° 50° 60°E 

Figure 2. Surface marine wind analysis from FNOC for I200Z August 6 1966. Isolmes repiesent contours of pressure at 4-mb intervals and define 
the location of high- and low-pressure systems. Add 1000 mb to all values < SO and odd 900 to all values > 50 to get correct pressure. H ind 
barbs are plotted every 5° of latitude and longitude in areas of change to show surface wind field. One barb = 1C kt (5 m/s) and one- half barb = 
5 kt (~ 2 m/s). 




SST as a “hot spot” when compared to surrounding 
ambient SST' values not impacted by the above condi- 
tions. As these favorable diurnal warming conditions 
continue, the mixed layer shallows and warms, creating 
a stratified layer that may destroy the surface acoustic 
duct. 

This so-called afternoon effect is critical to Navy 
ASW strategies (Urick, 1983, and Laevastu and 
Hubert, 1970). Knowledge of when the diurnal 
warming will negate use of the surface acoustic duct is 
important for a variety of sensor/weapon applications. 
The DOSL model could thus be used as one of several 
tools (along with the Thermodynamic Ocean Prediction 
System— TOPS) to help locate where the afternoon 
effect may take place. Appropriate action by naval 
forces could then lessen the chances that incorrect 
sensors or sensor deployment patterns were used within 
the specific domain of interest. 

Intense diurnal warming may create sharp, hori- 
zontal SST gradients that could be readily confused 
with other features related to real mesoscale ocean 
frontal boundaries. The generation of these gradients is 
especially disconcerting, since the regional FNOC three- 
dimensional thermal structure analyses (Optimum 
Thermal Interpolation System— OTIS) are gradually 
becoming more dependent on frontal “boguses” 
derived by analysts who extract frontal locations and 
strengths from infrared imagery. Thus, misleading SST 
gradient information could be passed onto the OTIS 
analysis if diurnal warming events are not screened out. 

Observations of sea surface warming due to 
insolation (diurnal heating) have been inadequate for 
many years (Cornillon and Stramma, 1985). Ship 
weather reports include SST readings that are often 
representative of the engine intake temperatures (Earle, 
1985). As such, they occur well below the water line 
(sometimes as deep as 15-20 m) and, thus, do not 
indicate the true SST. Some researchers have argued 
that this aspect is a positive factor and surface skin 
effects are negated, and that the ship values are a 
measure of the “mixed-layer” temperature. 

This logic does not hold true when dealing with the 
topic of diurnal surface heating. Surface heat fluxes 
and wind stress are the key to determining whether 
large temperature changes will occur at the ocean 
surface on a daily time scale. Therefore, the skin 
temperature (the apparent temperature of the sea that 
is in contact with the air) is the critical sea measure- 
ment, not some temperature value at a depth of 5 m 
or 10 m. This fact significantly raises the value of 
satellite infrared SSTs for the study discussed here. 

IV. Satellite-Derived Sea Surface 
Temperatures 

Satellite-derived SSTs have progressed significantly 
since the first infrared sensors were used. Meteoro- 
logical applications paved the way for the first series of 


environmental satellites, but oceanography was a 
distant by-product due in part to poor thermal 
sensitivity. This restriction was lifted when the 
TIROS-N series of National Oceanic and Atmospheric 
Administration (NOAA) polar oribiiers began using 
10-bit digitization (0.12°C sensitivity) over the whole 
range of possible SSTs. Small SST gradients could then 
be viewed and surface currents mapped as never before 
possible. However, the use of single-channel infrared 
data is limited by the attenuating effect of atmospheric 
water vapor. Maul (1983) clearly states the debilitating 
effect water vapor can have by calculating potential 
SST corrections ranging from nearly 0 to as much as 
10°C. These values are truly disconcerting when 
realizing that such errors can immediately rule out the 
vast majority of studies that rely on absolute SST 
values for quantitative studies. 

In 1981, NOAA began the operational production 
of MCSSTs by using the differential absorption 
observed in two infrared channels. The 10.3-1 1.3^m 
and 11.5-12.5/Km channels were used to generate a 
correction value applied to the 1 1 .5-12.5^m value that 
removed the major water vapor effects. Thus, the 
differing response between the two infrared spectral 
bands to the same columnar water vapor permitted the 
calculation of a correction factor, which fluctuates as 
the water vapor varies. 

The NOAA Advanced Very High Resolution 
Radiometer (AVHRR) is now a five-channel sensor 
that includes an additional spectral band between 
3.5fim and 3.93/im. This addition enables the use of 
a three-channel algorithm that can be superior to the 
two-channel, split- window equations (channels 4 and 5) 
when the ever-present noise in channel 3 is low just 
after launch or satellite outgassing. This problem will 
be corrected for the NOAA K,L : M series by 
incorporating slightly different sensor technology. 

The MCSST allows the adjustment to be done on 
a pixel-by-pixel basis and eliminates the need for 
additional data sources, such as atmospheric sounders, 
which are plagued with poor spatial resolution or 
radiosondes that provide sparse coverage over ocean 
domains. NOAA processes global area coverage 
(GAC) 4-km data in 2 x 2 unit arrays to produce 
8-km MCSSTs on a daily basis. The procedure to create 
MCSSTs has been refined over the years, and their 
accuracy has continually increased to the point where 
they now have an RMS error of 0.6-0.7°C (Strong and 
McClain, 1985) while producing 70,000 to 100,000 
global retrievals per day. This measurement milestone 
allows MCSST usage for a wide variety of research 
and operational projects. 

These National Environmental Satellite, Data, 
Information Service (NESDIS) MCSSTs are proc- 
essed orbit by orbit as the GAC data are downloaded 
from the spacecraft tape recorders. The time delay from 
sensor data collection to MCSST retrieval residing in 
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the NESDIS current day observation file varies from 
3 to 6 hours. These data can then be accessed by 
operational activities and incorporated into SST or 
three-dimensional ocean thermal analyses. 

This processing procedure is designed to create a 
reasonable global distribution to satisfy a variety of 
SST mapping requirements. However, not all cloud- 
free pixels are used to generate MCSSTs. The MCSST 
algorithm incorporates a search pattern that finds 
1 (low), 5 (medium) or 15 (high) observations per 
llxll GAC array (about 50 X 50 km). This factor 
limits the total number of MCSSTs produced and 
results mainly from several considerations — disk space, 
processing time, and overall data needs. Thus, 
hundreds of thousands of potential MCSST retrievals 
are not produced each day. 

The DOSL validation effort could use isolated 
matchups between day/night MCSST measurements 
and DOSL model output but is unrealistic for several 
reasons: no DOSL digital values were retained, only 
the contoured charts; interpolating coarse DOSL values 
to MCSST locations may induce large errors; and this 
method would not permit us to see whether the pattern 
of diurnal heating was being correctly reproduced. 
Therefore, MCSSTs (imagery) provides us with the 
synoptic scale coverage and highly accurate SST values 
required for DOSL comparisons. 


V. Validation Technique 

Cloud contamination severely restricts any validation 
effort that relies solely on infrared imagery. Therefore, 
instead of acquiring images in real time, the project 
relied on the DOSL model to indicate which regions 
should be viewed more closely. The first phase revolved 
around the collection of FNOC’s DOSL and marine 
surface wind charts for the 1986 Northern Hemisphere 
summer months. The charts were subjectively scanned 
to select only events that exhibited large diurnal 
warming over a spatial domain greater than 10° longi- 
tude (approximately 1000 km). 

DOSL charts for August, September, and most of 
October 1986 were scrutinized to find likely maximum 
diurnal heating events that would probably be 
associated with few clouds and thus indicate very good 
infrared imagery viewing conditions. Dozens of cases 
were located in the North Atlantic, the Mediterranean 
Sea, and the central and eastern North Pacific, as well 
as a few isolated occurrences in other spots. These 
potential case studies were then screened by taking into 
account the conditions usually found in each area. 

The North Pacific is well known for the extensive 
cloud coverage generated by the combined effects of 
the meteorological and oceanographic conditions that 
are manifested through air-sea interaction processes. 
These conditions are prevalent in the summer and the 
late fall when visible or infrared viewing is extremely 


poor at best. However, the DOSL model predicted 
large diurnal warming based on the forecast of light 
winds and little surface mixing even though skies were 
cloudy. The North Pacific was eliminated from further 
study, since the infrared data were not suitable. 

The Mediterranean Sea was initially thought to be 
a superior location for testing the DOSL model, but 
several problems arose. Excellent satellite viewing is 
the rule rather than the exception and permits regular 
infrared surveys. Light winds and insolation can quite 
frequently cause the model to forecast diurnal heating 
of 2°C. However, the regions of heating can be small, 
are quick to dissipate and reappear, and are difficult 
to time with regards to satellite overpasses. The region 
al'io frequently contains dust and aerosols, which can 
attenuate infrared signals and decrease the accuracy 
of measurements. 

Earlier efforts by Cornillon and Stramma (1985) and 
Stramma et al. (1986) showed that selected regions in 
the Sargasso Sea experienced diurnal warming in excess 
of I °C up to 30 percent of the time. DOSL charts were 
checked to locate days when significant heating was 
forecast and then nearly 100 NOAA AVHRR digital 
images were scrutinized to find cloud-free, day/night 
image pairs. These data were previously collected in 
real time by the NORDA Satellite Digital Receiving 
and Processing System (Hawkins et al., 1985) to 
support the GEOS AT Ocean Application Program 
(GOAP; Lybanon et al., 1987). GOAP combined the 
synoptic infrared coverage with precision GEOSAT 
altimeter sea surface height data to map mesoscale 
front and eddy features within the Gulf Stream system. 
All data were high picture resolution transmission 
data (1 km) from either Wallops Island, VA, or the 
NORDA tracking antenna. 

The DOSL model output indicated numerous cases 
where warming larger than l°C should have been 
visible in satellite imagery. However, a quick look at 
the digital images revealed very few suitable cloud-free 
pairs. Many day/night MCSST difference images were 
generated, but none were of sufficient quality for this 
study, even though cloud-free bits and pieces did verify 
well with DOSL values. 

Attention was then turned to the central and eastern 
North Atlantic, since the Azores-Bermuda high was 
often creating good cloud-free conditions. These 
conditions were evident in the persistently high DOSL 
diurnal heating forecasts. Numerous AVHRR digital 
images were ordered front NESDIS for further 
investigation. Both 1-km local area coverage and 4-km 
GAC data were included to see whether or not one was 
preferable to the other (note that Cornillon and 
Stramma, 1985, found 4 km sufficient for their study). 
Imagery for 16-18 and 28 August; 3, 15, 23, and 
30 September; and 2 and 7 October 1986 were then 
read in and initially scanned for cloud contamination. 
A quick survey concluded the August data were far 
superior, and they were then processed further. 
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The small number of cases were deemed insufficient 
in quantity, so attention was turned to other methods 
to increase the case study comparisons. Lessons learned 
during the 1986 summer were then incorporated into 
a near-real-time effort in the 1987 summer. DOSL 
charts were checked daily in the North Atlantic to try 
and catch the beginning of diurnal heating events. 
Conditions normally changed in logical phases as cloud 
systems with higher winds moved out and were 
followed by clearer skies and weaker winds. Direct 
communication between FNOC and NORDA 
personnel permitted the collection of several more 
cases. Each example, whether for 1986 or 1987, was 
processed in the following meaner: 

• NOAA-9 channels 4 and 5 data were calibrated 
and earth-located using software on NORDA’s Inter- 
active Digital Satellite Image Processing System 
(IDSIPS; Hawkins ct at., 1985). 

• Earth location was fine-tuned using European, 
African, and island landmarks and resulted in 
registration accuracies of 2 pixels or better. 

• Channels 4 and 5 were used for both day and night 
passes instead of 3,4,5 for the night and 4,5 for the 
day to be consistent and to reduce the introduction of 
noise inherent in channel 3. Conservative cloud tests 
were incorporated that likely erred on the side of 
eliminating some cloud-free pixels to make sure 
partially cloudy pixels were not allowed to slip through 
and introduce erroneous data points in the final results. 

• Multichannel split-window SSTs were calculated 
for all cloud-free pixels in both day and night imagery 
using the NESD1S operational equations valid for the 
period (Strong and McClain, 1985). 

• Daytime MCSST imagery was then subtracted 
from the nighttime versions to produce a day/night 
difference image. These imagery were then falsely 
colored in 0.5 P C increments to view more easily the 
regions of maximum diurnal heating. Clouds were 
masked in black or dark gray. 

The DOSL charts were not, however, in an easy 
format to compare directly with the images. Thus, the 
charts were digitized on IDSIPS and processed to 
the exact image map projection and scale. Isotherm 
contours at 0.2°, 0.5°, 1.0°, 1.5°, 2.0°, and 2.5°C were 
kept and all others were interactively erased to reduce 
the cluttered pattern in intense heating regions (DOSL 
output originally contains contours at 0.2°C intervals). 
DOSL contours could then be directly overlaid in a 
graphics plane on top of the falsely colored, day /night 
difference images. 


VI. Validation Results 

August 1986 was an excellent diurnal heating event 
time due to the meteorological conditions that existed 
in the central and eastern North Atlantic. The sub- 
tropical high dominated the central North Atlantic and 


contained vast areas of subsidence, as evidenced in 
pressure maps, such as shown in Figure 2. Infrared 
imagery was relatively cloud-free for a large region 
bounded by 25-40°N and I0-45°W. 

A. 16 Augusl 1986 

Figure 3 (I537Z 16 August 1986) is dominated by 
a warm, dark-red patch of SSTs oriented northeast to 
southwest off the extreme northwestern African coast. 
This section of very warm SSTs extends westward over 
20° longitude and is cloudless in the most intense area 
near 31 °N and 16° W. This very warm water is in direct 
contrast to the much cooler coastal upwelling along 
the northeastern African coast and contributes to a 
substantial SST gradient. 

The day/night MCSST image (Fig. 4) readily 
illustrates the scale of the diurnal heating taking place. 
A huge area of SST warming ranges from 0.0° to 0.5°C 
and covers a section surrounding the zone of maximum 
warming. A light-blue color was purposely chosen to 
offset it from the black and dark-gray shades used for 
cloud-contaminated pixels. This coloration permits the 
viewer easy visual access to the image content. (Note 
that there is little fading from one color to another in 
the color wedge, especially between 1° and 6°C, 
required to enhance each l.0°C day/night temperature 
warming bin). 

Close examination is needed to spot diurnal heating 
from 1 .0° to 2.0°C, even though it has been colored 
dark-blue to offset it from the previous light-blue color 
and the following dark-green (2.0° to 3.0°C) shade. 
This 1 °C bin is scattered about the image but is usually 
concentrated near the maximum heating sections. 

The 3.0° to 5.0°C diurnal heating bins (yellow and 
orange) are concentrated in limited sections and nicely 
outline regions of even higher heating. These values 
typically fall between 25° and 35°N and cover virtually 
the whole longitudinal range. This example of diurnal 
heating is impressive, even without the larger values 
found within the yellow- and orange-colored pixels. 
Note that conservative cloud tests were applied to 
ensure that few clouds contaminated the day/night 
images. Thus, some of the dark-gray pixels located near 
32°N and 23°W may actually be cloud-free or partially 
cloudy, but were masked out so that a simple test could 
be applied to the entire image. 

The DOSL contours (valid for the afternoon of the 
16th and overlain in white in Fig. 4) outline the general 
diurnal heating pattern. The very abrupt zone of 
heating just off the coast is well covered with a DOSL 
contour in excess of 1.0°C. The model forecast 
correctly predicts warming to occur around 25 °N 
across a wide ocean expanse. Two large pockets of 
2.5°C warming occur at 26°N, 25°W and at 32°N, 
35°W in the DOSL product, where the westernmost 
area agrees well with the day/night MCSST imagery. 
The eastern maximum is slightly misplaced to the 
southeast. 
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Figure 3. MCSST image at I537Z 16 August 1986 reveals extremely hot SST values off N W African coast. Clouds have been 
blacked out using a multichannel cloud-screening technique. SSTs have been falsely colored according to temperature scale 
in upper-right corner. Note patch in excess of 26 °C near 31°N, 16°W. 
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Figure 4. The 16 August 1986 day-night MCSST difference image, color coded according to temperature scale in 0.5°C and 
I.O°C increments. Light blue represents diurnal increases near l°C; yellow is near 3°C. White contours are FNOC DOSL 
SST change forecast values at 0.2°C intervals. Note how well elongated region off NW Africa correlates with DOSL local 
maximum. 
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The trend exhibited in Figure 4 represents excellent 
agreement between the DOSL output and the day/night 
MCSST imagery when considering the heating pattern. 
This agreement indicates that the wind stress and the 
surface heat fluxes are being correctly forecast by 
NOGAPS in terms of the overall meteorological 
conditions and are incorporated correctly by the DOSL 
model’s physics. The few regions of very intense 
heating, dramatically shown in imagery, reveal that the 
2.5°C upper bound within the DOSL model can be a 
limiting factor in certain cases. This threshold exists, 
since the model physics begin to break down for values 
above 2.5°C. 

One should also note the grid-point resolution 
involved in the input fields. NOGAPS has 2.5° spacing 
and is extremely mismatched when confronted with the 
realities of 1-km or 4-km imagery. Thus, the smoothing 
of heating effects is to be expected and will continue 
to be the norm until the atmospheric forcing and the 
DOSL model grid resolution are significantly reduced. 
In this light, the model appears to do quite well by 
mapping out the prime diurnal heating spots. 

B. 17 August 1986 

The following day, 17 August, also exhibited similar 
diurnal warming features, but not with quite the 
widespread intensity evident on August 16. Figure 5 
shows a large east-west section of DOSL-predicted 
afternoon warming that closely agrees with the 
0.0°-0.5°C (light blue) MCSST day/night difference. 
The eastern segment is not quite as sharply defined as 
in the previous day, but some cloud-free pixels have 
likely been masked out due to high surface reflectivity 
values caused by sunglint near 30°N, 15°W. 

The DOSL maximums near 30°N, 18°W and 33°N, 
32°W are somewhat cloudy; thus, good comparisons 
are difficult. However, the imagery clearly indicates 
several zones of 3° and 4°C heating between the two 
westernmost DOSL warming areas. Thus, the overall 
trend appears to be correct, but individual locations 
are hard to verify due to spotty cloud contamination. 

It should be noted that although the DOSL model 
has a 2.5°C threshold, the day/night MCSST 
maximum values of 3°-4°C compare well with the 
earlier results of Cornillon and Stramma (19i35) in 
the Sargasso Sea. Both areas are affected by the fringe 
of the Azores-Bermuda high, but on opposite sides. 
The eastern North Atlantic has the advantage that no 
major current system with large SST signals runs 
through the area. The Gulf Stream forces investigators 
to move to the south, away from the disturbing changes 
that can occur due to the energetic Gulf Stream cold- 
core rings. 

C. 6 August 1987 

This example involves a situation in the central North 
Atlantic using NORDA-acquired, 4-km GAC data. 


Figure 6 shows the data at about 4- kin resolution 
(considerably better than previous cases) and permits 
us to focus on a zone within the Azores-Bermuda high. 
A pronounced region of 0.0°-0.5°C warming is out- 
lined in light blue and agrees well with DOSL contours 
in shape and extent. More importantly, the main area 
of intense heating falls well within the 2.0°-2.5°C 
DOSL values. 

Diurnal heating values of 3°-4°C align along 
29°-30°N and 30°-45°W. These numbers exceed the 
DOSL threshold (as in the previous cases), but speak 
well for the model, considering the resolution for input 
forcing functions. The fact that clouds lie to either side 
of the maximum warming indicates that NOGAPS 
appears to have done well in forecasting the 
meteorological conditions. The extra region of high 
values in Figure 6 suggests that NOGAPS is in error. 
A band of clouds oriented northeast-southwest 
separates the two high-value sections, but DOSL misses 
the northwestern segment by labeling it as fi,2°-0.5°C 
instead of the 2°-4°C heating present in the imagery. 

The corresponding surface winds at 00Z 6 August 
1987 are outlined in the FNOC Marine Wind Analysis 
(Fig. 7). This analysis is separate from the NOGAPS 
fields and includes a combination of ship-, island-, 
and satellite-derived wind retrievals. Isobars at 
4-millibar (mb) contour intervals are superimposed 
with high- (H) and low- (L) pressure centers and with 
vectors (half barb = 5 kt, full barb = 10 kt, etc.). 

As expected, the Bermuda High dominates the 
meteorological pattern throughout most of the North 
Atlantic from 25°N-35°N. Winds near 15-20 kt 
surround the region, but much calmer conditions reside 
in the high-pressure center as warm, dry air subsides 
from above. This situation helps to reduce tropical 
cloud development and promotes diurnal heating. 

These light winds are also verified by the Special 
Sensor Microwave/Imager (SSM/I) launched June 
1987 aboard a Defense Meteorological Satellite 
Program spacecraft. This passive microwave sensor is 
similar in frequency and purpose to the Nimbus-7 
scanning multichannel microwave radiometer, but has 
finer spatial resolution (25 km) and better on-board 
calibration. It can measure surface wind speed (2 m/s), 
columnar water vapor (0.3 g/cm 2 ), and sea-ice 
concentration (10 percent), among other parameters 
(Hollinger, 1989). 

The SSM/I can measure surface winds in essentially 
all weather conditions. Heavy rain cells can cause 
erroneous retrievals, but these values can be detected 
by flagging regions of heavy precipitation, another 
parameter sensed by the SSM/I. The 1394-km SSM/I 
swath permits wide sector viewing during its 14 polar 
orbits. However, data gaps exist below 50°, since 
consecutive passes do not overlap. The gap increases 
to 1000 km at the equator and is then filled in largely 
by orbital data 12 hours later. Complete global 
coverage occurs every 36 hours. 
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Figure 5. The 1 7 Augu i 1986 day-night MCSST difference image , color coded according to temperature scale. General agreement 
in western portion , but DOSL values are too high in area of previous maximum , as seen in Fig. 4. 
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Figure 6. The 6 August 1987 day-night MCSST difference image , color coded according to temperature scale. Iris image focuses 
on area near center of subtropical high to reveal intense diurnal heating in excess of 4°C, but is well mapped by the DOSL 
contours. Note this figure uses higher spatial resolution than incorporated in Figs. 4 and 5. 
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LATITUDE 



Figure 7. NOAA National Hurricane Center surface weather analysis for I200Z 16 August 1986. Pressure contours at 4-mb 
intervals; wind vectors of one- half barb = 5 kt (- 2 m/s) and one barb = 10 kt (5 m/s). The two sets of numbers next to 
each ship weather observation are air temperature (°C) and surface pressure ( multiply value by 0. 1 and add 1000 mb to get 
correct value). Sky cover is marked by the black portion of observation circle (i.e., open circle with a tine denotes no clouds, 
half a black circle indicates presence of 50% clear skies). Note dear skies in center of subtropical high near 30° N, 40° W. 



LONGITUDE 


Figure 8. SSM/I surface winds on 6 August 1987 for 0700Z to 
IOOOZ. The 25-km wind retrievals have been mapped into a two- 
dimensional field to examine the general wind speed patterns. 
A large zone of winds less than 10 kt (5 m/s) parallels 30°N. 


Figure 8 represents several passes covering the DOSL 
image domain for 5 August 1987 (no SSM/1 values 
available on 6 August due to satellite coverage gaps). 
The wind speed retrievals are contoured at 5-m/s 
intervals and are shaded to help the viewer easily 
recognize the zones of constant wind speed. Light 
winds dominate the region and provide accurate 
backup information to the marine wind analysis. 

VII. Future Considerations 

The NOGAPS model has undergone considerable 
revision since 1987 and now contains an entirely 
new treatment for the atmospheric boundary layer. 
The new approach has proven to be far superior to 
the earlier version, so heat fluxes for DOSL input have 
been upgraded significantly since this study began. 

The NOGAPS and DOSL spatial resolution of 
2.5° is quite coarse when attempting to nail down 
diurnal heating features that can be very localized. 
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Future reductions in model grid spacing will help to 
solve this problem. 

Sampling day and night SSTs to validate DOSL 
results with one polar orbiter is difficult. Any follow-on 
studies should seriously consider the use of 
geostationary data that provide views every 30 minutes 
and that remedy some of the cloud contamination 
problems caused by compositing over short intervals 
as clouds move. The critical element is the accuracy 
of SSTs from the Geostationary Operational 
Environmental Satellites (GOES). The V1SSR 
Atmospheric Sounder has reached accuracies of less 
than 1°C, but improvements are needed before it is 
equivalent to MCSSTs. This will likely come about with 
GOES-NEXT in 1990-91, where multi-infrared 
channels and 10-bit digitization will open up more 
applications for these platforms. 


VIII. Summary and Conclusions 

As shown by the work of Cornillon and Stramma 
(1985), Stramma et al. (1986), and many others, diurnal 
SST variability is ubiquitous and widespread during 
the spring and summer. Forced by day/night variations 
in the net surface heat flux and strongly modulated 
by the magnitude of the surface wind stress, this 
phenomenon dominates short-term variability in upper- 
ocean thermal stratification (Price et al., 1986), and 
thus critically impacts on surface-duct acoustics. 

Diurnal SST variability is also important when using 
and interpreting satellite SST data. Satellites measure 
the skin temperature of the ocean, which may or may 
not be characteristic of the average temperature above 
the seasonal thermocline, depending on the amplitude 
and phase of the diurnal SST cycle at the time of 
measurement. Furthermore, because the diurnal SST 
response is such a strongly nonlinear function of 
the wind speed, smoothly varying patterns of the 
synoptic wind field tend to produce sharp horizontal 
gradients in diurnal SST features, as demonstrated by 
the satellite imagery presented here. These sharp 
horizontal gradients would be easy to misinterpret as 
water mass bou idaries. 

Price et al. (1986) developed scaling relationships 
relating the amplitude of the upper ocean’s diurnal 
response to wind-stress and heat-flux parameters. They 
verified these relationships, referred to as the DOSL 
formulation, by using observed winds, heat fluxes, and 
SSTs at a moored buoy. 

Results from a real-time synoptic implementation 
of the DOSL model at FNOC were presented and the 
predictions of this model were verified with satellite 
MCSST data. The case studies discussed here were 
biased to the extremely warm diurnal heating events, 
since the infrared imagery must obviously be cloud- 
free for use and, thus, large-magnitude diurnal heating 
has an increased chance to occur. The DOSL correctly 


defined the shape or pattern of each heated region, but 
the upper bound on the warming is too low by virtue 
of the self-imposed upper bound of 2.5°C. Day/night 
MCSST differences reached 3°-5°C in several cases 
when light winds combined with intense summertime 
insolation in the subtropical highs. 


IX. Recommendation 

DOSL has proven to handle these extreme cases well 
in a qualitative manner and thus can be used for several 
immediate applications. This report recommends that 
the DOSL model be declared operational and used to 
fill the gap existing in forecasting daily SST diurnal 
warming events. 
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